^ , 


s 
s 

I 

O 

o 


> 
in 
oo 

00 

o 

0^ 
OS 
-I— > 


I 

O 


X 


Comment on "Density Functional Simulation of a 
Breaking Nanowire" 

In a recent Letter [Q, Nakamura et al. described first 
principles calculations for a breaking Na nanocontact. 
Their system consists of a periodic one-dimensional array 
of supercells, each of which contains 39 Na atoms, orig- 
inally forming a straight, crystalline wire with a length 
of 6 atoms. The system is elongated by increasing the 
length of the unit cell. At each step, the atomic con- 
figuration is relaxed to a new local equilibrium, and the 
tensile force is evaluated from the change of the total 
energy with elongation. Aside from a discontinuity of 
the force occuring at the transition from a crytalline to 
an amorphous configuration during the early stages of 
elongation, they were unable to identify any simple cor- 
relations between the force and the number of electronic 
modes transmitted through the contact. An important 
question is whether their model is realistic, i.e., whether 
it can be compared to experimental results |Q obtained 
for a single nanocontact between two macroscopic pieces 
of metal. In this Comment, we demonstrate that with 
such a small unit cell, the interference effects between 
neighboring contacts are of the same size as the force 
oscillations in a single nanocontact. 

In order to understand how the close proximity of the 
nanocontacts in the model of Ref. [|l| may alter the ener- 
getics of the system, we consider a system of two identical 
nanocontacts in series, connecting two macroscopic wires. 
We model the metallic nanocontacts as constrictions in 
a free electron gas, with hard- wall boundary conditions, 
and obtain the energetics of the system from the elec- 
tronic scattering matrix The scattering matrix of 
the compound system may be obtained as a geometric 
series in the scattering matrices of the individual con- 
tacts (which are taken to be symmetric under inversion, 
for simplicity), while the scattering matrix of a single 
contact may be evaluated using the adiabatic and WKB 
approximations [p|, which are quite accurate for contacts 
of smooth shape . The total grand canonical potential 
of the system is found to be the sum of the contributions 
of the individual contacts, plus an interference term 


An 


dEf{E) ^ tan" 


R,{E) sin[2e,{E)] 
l + R,,{E) cos[2e^,{E)y 


where f{E) is the Fermi-Dirac distribution function, and 
Ru{E) and 0^{E) are the reflection probability and scat- 
tering phase shift, respectively, of the i^th electronic mode 
for a single nanocontact. 

The magnitude of the correction to the cohesive force 
in the supercell arrangement of Ref. jl] arising from in- 
terference effects between neighboring supercells is AF — 
— 9[Ar2]/9Lccii, where LcoU is the unit cell length. Inter- 
ference between more widely separated supercells would 
lead to an additional correction. Fig. 1(b) shows that 
for the unit cell size considered in Ref. |Q| (Lccii = 17 — 
31 A = 2.5 — 4.5Ai;'), the interference correction to the 


cohesive force is comparable to the force oscillations of 
an individual nanocontact. For comparison, the conduc- 
tance of a single nanocontact and the interference correc- 
tion thereof are shown in Fig. 1(a). For a single contact, 
there is a clear correlation between the conductance steps 
and the force oscillations. However, the large interference 
correction would strongly suppress any correlations be- 
tween the force calculated in the supercell arrangement 
of Ref. and the conductance of a single contact. 

In order to explain the correlations between cohe- 
sion and conductance observed experimentally in metal- 
lic nanocontacts 1^ , it is essential to treat the energetics 
and transport of the system on an equal footing. This has 
been achieved in our free-electron model Q . The interfer- 
ence term scales as AF C(-^cx!ii) [since 9v{E) oc iccii], 
so it would be worthwhile to perform larger-scale "first 
principles" simulations to address this question. 
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FIG. 1. Conductance and force for a metallic nanocon- 
tact with dimensions comparable to that of Ref. [1]: initial 
radius=AF, initial length=2.5AF. (a) Conductance for a sin- 
gle contact (solid curve), and interference term (dotted curve, 
offset by 3); (b) force oscillations for a single nanocontact [3] 
(solid curve), and interference term AF = — 9[Ar2]/9Lccii 
(dotted curve). 


